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A Blume-Emery-Grifiiths perceptron model is introduced and its optimal capacity is calculated 
within the replica-symmetric Gardner approach, as a function of the pattern activity and the imbed- 
ding stability parameter. The stability of the replica-symmetric approximation is studied via the 
analogue of the Almeida-Thouless line. A comparison is made with other three-state perceptrons. 



PACS numbers: 02.50-r, 64.60.Cn, 75.10.Hk, 87.18. Sn 



I. INTRODUCTION. 

Recently an optimal Hamiltonian for a multistate net- 
work has been put forward [|| on the basis of in- 
formation theory by maximizing the mutual information 
content of the system. For a two-state network, this 
Hamiltonian equals the well-known Hopficld Hamiltonian 
extensively studied in the literature || , || . For a three- 
state network one finds a Blume-Emcry-Griffiths (BEG) 
spin-glass type Hamiltonian ||. As spin-glasses these 
models have been studied for some time now. Thermo- 
dynamic as well as dynamic properties are discussed in 
the literature for disorder in both the quadratic and bi- 
uadratic interaction. Many references can be found in 
As a neural network model its study has been started 
only recently Q , |I(J . But it turns out already that both 
the maximal capacity and the basin of attraction of this 
network are enlarged, at least for Hebb rule learning, in 
comparison with the standard three-state networks like, 
e.g., the Q-Ising spin-glass models. 

A natural question is then whether these improved re- 
trieval quality aspects are restricted to the use of the 
Hebb rule or whether they are intrinsic properties of the 
model. In the same context, a further question is then 
whether we can extract a perceptron type model with an 
optimal performance out of this BEG recurrent network. 
The perceptron is by now a well-known and standard 
model in theoretical studies and practical applications 
in connection with learning and generalization H, [Q, 
- p. Consequently, a number of extensions includ- 
ing many-state, graded response and colored perceptrons 
have been formulated in the literature |Ij]-p8|. 

The aim of this work is precisely to introduce such 
a BEG-perccptron model and, in particular, to study 
its Gardner optimal capacity. Although the method for 
doing that is standard and well-know by now [l9|] , ]2C[] 
its generalization to the problem at hand is highly non- 
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trivial. Nevertheless we have succeeded in obtaining a 
closed expression for the replica symmetric approxima- 
tion to the Gardner optimal capacity. 

The paper is organized as follows. In section 2 we 
recall the BEG Hamiltonian and define the BEG percep- 
tron model. Section 3 presents a closed analytic formula 
for the replica-symmetric Gardner capacity of this model 
and studies its behaviour as a function of the imbed- 
ding constant and the activity. Comparisons with other 
three-state perceptrons are made. In section 4 the sta- 
bility of the replica symetric solution is studied using an 
extension of the de Almeida-Thouless analysis. The an- 
alytic form of the two replicon eigenvalues is obtained. 
Stability is found to be broken for smaller values of the 
activity and for very small imbedding stabilities. Section 
5 presents some concluding remarks. In the appendices 
further technical explanations are given. 



II. THE BEG PERCEPTRON. 

Consider a neural network consisting of TV neurons 
which can take values <7j, i — 1, . . . , N from the discrete 
set S = {—1, 0, +1}. The p patterns to be stored in this 
network are supposed to be a collection of independent 
and identically distributed random variables (i.i.d.r.v.), 
I , . . . , p with a probability distribution 



PiO = X - 1) + n*(tf + 1) + (i - WW) (i) 
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with a the activity of the patterns so that 
1 
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N^oc 
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Given the network configuration at time t, <jn = 
{cj(t)},j — 1, . . . , AT, the following dynamics is consid- 
ered. The configuration <tat(0) is chosen as input. The 
neurons are updated according to the stochastic parallel 
spin-flip dynamics defined by the transition probabilities 



Pr (cn(t + 1) = s' eS\<r N (t)) = 



exp[-(3ei(s'\cr N (t))} 

E s es ex P[-^ e i( s l CT ivO))] 
(3) 
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Here the energy potential ej[s|<rjv(i)] is denned by 

ei[s\(T N (t)] = -shi((T N (t)) - s 2 9 l {a N {t)) , (4) 

where the following local fields in neuron i carry all the 
information 

h N ,i(t) = Jiyo-j(*). = K ^°)^) ( 5 ) 

with the obvious shorthand notation for the local fields. 
For synaptic couplings Jy and of the Hebb-type 



J%j — 
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(6) 



E((£f) 2 -a)((^) 2 -a) (7) 



a 2 (l-a) 2 iV ^ v 
the corresponding neural network Hamiltonian 

fl" = ~\ J] Jy CT^ - - X Xyff?Cr| , 



(8) 



has been discussed recently [||. It has been found that 
the capacity and basin of attraction has been enlarged in 
comparison with other three-state networks. 

We would like to understand whether these better 
retrieval quality is an intrinsic property of the model. 
Therefore, we want to answer the following question: 
given the set of p patterns specified above, is there a net- 
work (the best possible network of the BEG-type) which 
has these patterns as fixed points of the deterministic 
form of the dynamics considered above? At zero tem- 
perature the updating rule of this dynamics @-([|) is 
equivalent to the gain function 

ai(t + l) = sign(h N ,i{t))@(\hN,i(t)\ + 9 Nti (t)) 

= g(h N>i (t),e N>i (t)) (9) 



and {J, K} = { Jj, Kj} denoting the configurations in the 
space of interactions. The factor N' 1 / 2 is introduced to 
have the weights Jj and Kj of order unity. 

The aim is then to determine the maximal number of 
patterns, p 7 that can be stored in the perceptron, in other 
words to find the maximal value of the loading a = p/N 
for which couplings satisfying (ji^|)-( |Tl|) can still be found. 
Following a Gardner-type analysis ]19| the fundamental 
quantity that we want to calculate is then the volume 
fraction of weight space given by 

f p 
V = J dJdKp(J,K) JJxeff(^^;«) (12) 
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with the characteristic function 
X${h' J ',6 fi ;K) =5^^^) 

= (^) 2 e(|^| + r- K )e(^ 
+(i-(^) 2 )©(-l^l-^-«) 



(13) 



where n is the imbedding stability parameter measuring 
the size of the basin of attraction for the ^-th pattern and 
p(J,K) is the following normalization factor assuming 
spherical constraints for the couplings 



p(J,K) 



5(3-3- N)S(K -K-N) 
J x x d3dKS{3 ■ J - N)6(K -K-N) 
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In order to perform the average over the disorder in the 
input patterns and the corresponding output we employ 
the replica technique to evaluate the entropy per site 



lim — ((lnVO) 



(15) 



where ((•■•)) denotes an average over the statistics of 
inputs and outputs recalling (III). 



with the Heaviside function. Considering the percep- 
tron architecture (N inputs with couplings Jj and Kj and 



1 output) we say that a given pattern, = 1, 
is stored if there exists a corresponding output 



with 
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III. REPLICA SYMMETRIC ANALYSIS 

In the replica approach the entropy per site v is com- 
puted via the expression 

v= lim lim (((V 71 } - l) = lim lim— In (V n } 

(16) 

(11) where V n is the n-times replicated fractional volume 



.,N, 
(10) 



, n II n p 

J [i[d3 a dK a 5(3 a ■ 3 a - N)s(K a -K a - n)}(( ]J 



(17) 
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whereby we can forget, since the couplings are contin- 
uous, about constant terms such as the denominator in 
( |l4| ) . The calculation then proceeds in a standard way al- 
though the technical details are much more complicated. 
For a short account we refer to Appendix A. Here we re- 
strict ourselves to the following important remarks. The 
main order parameters appearing in the calculation are 



<? Q /3 = ^J Q -J /3 , r af3 = -J-K Q • K' 3 , a</3 
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L a = -pViT" Va. 



(18) 



Of course, in the replica symmetric (RS) approximation 
we are focussing upon here, q a p — q, r a p = r,L a = L. 
The first two order parameters are the overlaps between 
two distinct replicas for the couplings J and K, the third 
one arises from the fact that the dynamics (j^) and, hence, 
also the characteristic function (13), contains a second 



field 9, quadratic in the patterns. We remark that it 
describes the relative importance of the active versus the 
non-active neurons. Actually, in the calculation ah will 
be the important quantity with a the second moment of 
the pattern distribution, i.e., the pattern activity. 

The RS optimal Gardner capacity is obtained when 
the overlap order parameters q and r go to 1. It is clear 
that these limits have to be taken simultaneously but, 
in general, their rate of convergence could be different. 
Therefore, we introduce (1— r) = 7(1— q) where 7 is a new 
parameter which one also needs to extremize. We expect 
this parameter 7 to depend on the pattern distribution 
through the activity a. 

Pursuing this approach then leads to 



ajisia, k) = — extr lim 



1 + l/ 7 



T^'q^i 2(1- q)gf s (q, 1 ,L) 
where gf 5 (g, 7, L) reads 



(19) 



g? s (q, 1 ,L)= I V(h )V (y^8o - I) (( In 



dh 



dO 



O s y/2n(l - q) y/2n(l - q) 
I 



exp 



(h - ho) 2 + (6 - 6 f 



2(1-9) 



(20) 



with I = aLj i/op — a), where T>(ax + 6) = 
(27r) _1 / 2 adxexp[(— 1/2) (ax + b) 2 ) and where the in- 
tegration region fl^ is determined by the Heavi- 
side functions appearing in the characteristic function 
X£(i/a/i, ^/7a(l — a)9; k) defined in (|D|). The expression 



( pp| ) for the function g± s suggests that an asymptotic ex- 
pansion to compute the limit q — ► 1 is possible. Indeed, 
after some tedious algebra (see Appendix B) wc find for 
this limit 



9? S (q,J,L) 



3 r 

' ~20^qj S j n V ( ko + «/V5©(V7«o - <Cl n (h 0> 6 ) 
^TT) ? J n , V (WVl0o ~ «) d K ^0, Bo) + 0(1/(1 q)) 



(21) 



with u = (aL + n)/ \J a(l — a). The integration regions 
read 



Tit 

n 2 
K 

K' 2 



and the corresponding integrands are given by 
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FIG. 1: The optimal capacity ors as a function of the pattern 
activity a for several values of the stability constant k. The 
dots at a = 1 refer to the optimal capacity of the two-state 
perceptron. 



FIG. 2: The optimal capacity curs as a function of the stabil- 
ity k for several values of the pattern activity a. The straight- 
dotted line corresponds to the optimal capacity of the two- 
state perpectron. 



with 7' = \/7(l — a) and where we remark that the d m i n 
are minimal distances between a point in the different 
integration regions lZi,lZ'i,i — 1,2,3 and the border of 
f^(see Appendix B). This may allow for a possible geo- 
metrical interpretation of the Gardner optimal capacity 
in the space of local fields as it has been suggested for 
the Q-state clock model in pl[ . 

After inserting (^l|)-(^3|) in ( |l9| ) and extremizing nu- 
merically with respect to L and 7, we find the results 
presented in figures In fig. 1 the capacity aiis 

versus the activity a is shown for several values of the 
imbedding stability constant k. For bigger k, the capac- 
ity becomes, of course, smaller. For a = 1, i.e., binary 
patterns, we find back the original Gardner results, as we 
do in fig. 2 showing aus as a function of k for several 
values of a. Smaller activity indicating a growing pres- 
ence of zero-state neurons leads to bigger capacities. Of 
course, this does not mean a priori that also the informa- 
tion content of the system is increased. For completeness, 
we remark that the parameters I = ahj \f a{l — a) and 7 
that we have extremized over, depend rather strongly but 
smoothly on the pattern activity. For a = 1 we find back 
the two-state perceptron value for L, i.e. L = 0(7 = 00), 
and 7 = 00. Finally, in order to have an idea about the 
information stored into the network we plot in fig. || the 
information content per neuron 



7(a) = - 



a RS 

In 3 



aln(-) + (1 - a)ln(l - a) 



(34) 



For a = 1 our result is again consistent with the simple 
perceptron result |l9| . Comparing with other three-state 
neuron perceptron models we recall that for n = and 
uniform patterns the Q = 3 Ising perceptron can maxi- 
mally reach an optimal capacity equal to 1.5, depending 
on the separation between the plateaus of the gain func- 
tion (see JL4|, M) for the precise details) and the Q = 3 




FIG. 3: The information content per neuron, /, as a function 
of a, for k = 0; 0.2; 0.4; 0.6; 0.8 (from top to bottom) 



clock and Potts model both reach an optimal capacity 
of 2.40 while the value for the BEG perceptron 

found here is 2.24. Here we have to recall that the Q = 3 
Ising perceptron and the BEG perceptron have the same 
topology structure in the neurons, whereas the Q = 3 
clock and Potts models have a different topology. 



IV. STABILITY OF THE REPLICA 
SYMMETRIC SOLUTION 

From the work of Gardner |l9| we know that for the 
binary neuron perceptron the RS solution is marginally 
stable against RS breaking (RSB) fluctuations. From the 
work on multi-state Q-Ising neurons pi[ we know that 



■5 



the RS solution may be stable or unstable depending on 
the gain parameter, the number of spin states and the 
distribution of the patterns. Furthermore, in general, 
increasing the imbedding stability parameter k lowers the 
capacity and enhances the stability against RSB. Using 
these results for the Q — 3 spin states as a guide we 
also expect breaking for the BEG perceptron model at 
hand. To confirm this and find out the precise interval 
of a values where breaking occurs^ we generalize the de 
Almcida-Thouless analysis 22 



First, the hessian matrix associated with the func- 



tion eq. (A8), is computed, and then the eigenvalues 
are determined. As usual, two types of eigenvalues are 
found: longitudinal eigenvalues describing fluctuations 
within RS and transverse eigenvalues describing stability 
against RSB. We find four transversal eigenvalues each 
with degeneracy ^n(n — 3). In the limit q — > 1 they can 
be calculated explicitly in terms of the minimal distances 
occuring in (p8|)-(|33|). The result reads (for more details 
we refer to Appendix C) 



A+ = ~(A, + A r ) + y(A q A r f + 4A2 (35) 
A_ = \(A q + A,.) - \^(A q - A r f + 4A* (36) 

T+ = 2(A2 - A A ) i Aq + Ar + (A? + A?)(A " ~ A?Ar) + \/ 4A ' + [ A « - A r + ( A § - AO(A,A r - A?)] 2 } (37) 



2(A 2 

X A A ) {A, + A r + (A ? + A ? )(A 2 - A,Ar) - ^Al + [A, - A,. + (A 5 - A f )(A,A r - A?)] 2 } (38) 



T- = 



2(A 2 



c ^q^r) 



with the A's given by 



a<? = ^^E J n ^ + ^/^^-t){\^d^ m {hM\ 



(1 (1 ^ g "f E j n , V ^ V7^o - «){^ J^™L(V )} 2 + o(l/(l q)) (39) 



A r = 



TZ^f E f K V ^ V7fo - "){\^$L<te> e,)) 2 + o(l/(l q)) 
± £ V(H + - ,){ | ^^(fc, , )} 2 

£ V(h , V7^o - u){ g Qf^ffiin^ ^} + ^/(l - «)) 



7 2( 

(1 - a)a RS 
' 7 2( 



(1 - a)a RS 



(40) 



(41) 



7(1 -9) 2 

A ? =(1-<Z) 2 (42) 
A ? = {l-rf =1 2 {l-qf . (43) 



Then two replicon eigenvalues, X Rl and \r 2 , can be de- 
fined as 

A fll = A+r_ A fe = A_r+ (44) 
Stability of the RS solution requires that both X Rl , A _r 2 < 



0. In fig. |j-|6| we present the numerical results concern- 
ing the stability analysis. In fig. |] the first replicon 
eigenvalue \ Rl is shown as a function of a for several 
values of k. It is seen that for small values of k this 
eigenvalue becomes positive for smaller values of a and 
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FIG. 4: The first replicon eigenvalue Arj as a function of a 
for several values of k. The dots at a = 1 refer to the optimal 
capacity of the two-state perceptron. 



hence replica symmetry is broken. We remark that for 
a = 1 our results are consistent with those of Gardner 
p9| . Fig. U presents a closer view of this for k = 0. For 
< a < 0.48(8) the RS solution is unstable. Storing 
only zero-state spins, a = 0, or binary spins a = 1 leads 
to marginal stability. As a first explanation one could 
remark that for increasing a, allowing more ± states, the 
disorder is increased up to about a uniform distribution 
of patterns, a — 2/3. It is clear that for bigger k, the 
stability against RSB increases. In fact for k > 0.0061 
already no more breaking occurs. Finally, fig. 6 shows 

0.01 i — — . — — . — — . — — . — — i 



0.00 
fl, -0.01 
-0.02 




-0.03 1 — ^ — 1 — ^ — 1 — ^ — 1 — ^ — 1 — ^ — 1 
0.2 0.4 0.6 0.8 1 

a 

FIG. 5: The first replicon eigenvalue as a function of a 
for k — on a different scale. RSB occurs for smaller values 
of a. 

that Xr 2 is always negative and, hence, plays no role in 
the breaking of the RS stability. 



FIG. 6: The second replicon eigenvalue \r 2 as a function of 
a for several values of k. 

V. CONCLUDING REMARKS. 

In this work we have introduced a perceptron model 
based upon the recently studied Blume-Emery-Griffiths 
neural network, containing ternary neurons. We have 
obtained an analytic formula for the replica symmetric 
optimal Gardner capacity. For the imbedding stability 
constant equal to zero and uniform patterns, e.g., we 
find a bigger optimal capacity, curs = 2.24, than the one 
for the Q — 3 Ising perceptron, a us — 1.5, which has 
the same topology structure for the neurons. Since, in 
general, perceptrons turn out to be very useful models 
in connection with learning and generalization this is an 
interesting observation. It is also consistent with earlier 
results derived for the Hebb rule. 

We have studied the stability of the replica-symmetric 
solution by generalizing the de Almeida-Thouless analy- 
sis and deriving an analytic expression for the two repli- 
con eigenvalues that play a role in the Gardner limit. 
Breaking only occurs for small activities and very small 
imbedding constants, k < 0.0061. This is consistent with 
the stability results found for the Q — 3 Ising percep- 
trons. 

These results strenghten the idea that the better re- 
trieval properties found for the Blume-Emery-Griffiths 
model in comparison with the Q = 3 Ising model are not 
restricted to the specific Hebb rule but are intrinsic to 
the model. 
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APPENDIX A: REPLICA ANALYSIS AND 
REPLICA SYMMETRIC ANSATZ. 

In this appendix we outline the main steps in the cal- 
culation of the n-times replicated volume (ujb extending 



J 



|19| to the case at hand. In order to perform the quenched 
average we use the (^-function representation 



00 dhf,dh" 



2tt 



00 deye? 



exp 



3=1 

, N 



to take the local fields out of the characteristic function and obtain 

\\ q=1 fi=l II J a=l ft=l 

II n p 

(nn«p 

\\ a =l fi=l 



V\2 



n p 



[ II II ^r^r] ex p [* E E + ^ 



M\2 



a — 1 /1— 1 

\\ a=l ^=1 // £g 

Introducing the order parameters ([l8|)and their conjugate variables, and using the identities 

dq a pdq a p 



■ft* n g a N 



1 = 



n 

1 a</3 



2m /N 



■ exp 



q a p[Nq a p- J" • J' 5 



n dr aj3 dr al 3 
2m IN 

' a<(3 1 



exp 



7Vr Q/3 — K a ■ K' 3 



1 



n 



dL a dL a 



exp 



* 27r/ViV 

allows us to express the replicated fractional volume as an integral over them, viz. 

' x ' ' dL a dL a ^ 1 7 f dE a dF a 1 |" -p-j- dq af} dq af) dr af} dr afi 

with <E> given by 



(V) 



oc 



r -pr efi/*eUv"-| r -j-j- dii" &f<'"-| r -pr 



a</3 



exp 



$ — aG\(q a p, r a p, L a ) + G2(Q a f3, R a p, L a ) + G3(q a p, r a p,Q a p, R a p) 



where 



d = In 



n 

[n 



2tt 2tt 



n n n 

exp [i £ (/i Q /i a + 6 a e a ) - ia 0° L " " - E 



a,/3=l 



a(l — a) 



G 2 



G, = 



hi 



j J| dJ a dK a 



exp 



Q = l 

1 

~ 2 



a,/3=l 



1 " 



a, /3=1 



r 



(Al) 

(A2) 



(A3) 

(A4) 
(A5) 

(A6) 

(A7) 
(A8) 



(A9) 
(A10) 
(All) 



and 



Q tt/ 3 = £ a <W + g a /j(l - <W) (A12) 

i? Q /3 = F a $af3 + f a/3 (l " « a/J ) (A13) 
Qa/3 = <W? + - <W) ( Al4 ) 



We remark that the 5-function representation of the local 
fields has allowed us to perform the calculations until this 
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point without using an explicit form for the characteris- 
tic function x^(ft. Q , 9 a ; k). Using the RS ansatz 4> can 
be simplified further and the saddle-point equations for 
Q,Q,R,R become algebraic so that they can be solved 
explicitly, leading to the result (119)- (EOj) . 



APPENDIX B: q -> 1 LIMIT 

In order to compute the asymptotic expansion of ( pQ ) 
we proceed as follows. We split the integral over (ho, 9o) 
into two parts, i.e., determined by the Heaviside func- 
tion in X(i an( i hs complement C(Q^). The first integral 
gives zero contribution in the limit q — > 1, while the sec- 
ond one gives a contribution of order (1 — q)~ x . Indeed, 
the integration over (h, 9) parametrized by q is nothing 
but an exponential Dirac-delta representation. When- 
ever the peak of this delta representation lies in the re- 
gion f2£, which means that (ho,9o) € fi^ the integral 
results in a finite contribution. The contributions of or- 
der (1 — q)^ 1 arises from the points (ho,9o) £ C(Sl^). 
Therefore, we can write 



y ' 






Tlx 






(h ,9o) 










^ ► 

X 


(h ,e ) / 




K 3 



(h ,e ) 

FIG. 7: Schematic representation of the subregions and min- 
imal distances for C(Q^ = 1) 



V(h )V(^9 -l)Hl} s (h ,9 Q ) 



where we have introduced the shorthand notation 
f dh d9 

' Jn, y/2n(l-q) 727r(l-g) 



(Bl) 



mdho,9 ) 



x exp 



2(1 -q) 



.(B2) 



Next, for a given (Jiq,9q) £ C(Q^) the main contribution 
arising from the function [1]^(/iq, 9q) is obtained for those 
points (h, 9) S which minimize the distance (h— ho) 2 + 
(9 — 9q) 2 . To calculate this minimal distance, we split up 
C(r^) into three subregions according to fig. ^in the case 
of £ = 1 



f ^0 < -7= 

1 \e > o 



(B3) 

n 2 = {( h °-7z)^ r ^ <9 ° <0 (B4) 
\ h o<t« 
9 <0 



n a = 



-4= < ho < -7= 



9oV7(l- - (B5) 



Computing the minimal distances for such subregions 
is straightforward and leads to 




{ho, Oq 



(h ,e ) 



FIG. 8: As figure @ for C(Q e = 0) 



d Ul = 

min 



jU 2 



h 



(B6) 
(B7) 



l + ,(l-a)KV = a- e ° VW ^ r) - h0 ) (B8) 



By redefining h a — n/^fa — ► ho and 7' = 1/7(1 — a) we 
recover the expressions (p8|)-(|3(i|). 

We proceed analogously for the region C(fl^ = o). We split 
this region into three subregions as shown in fig. 
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K 
K 
K 



h >0 

~^=—,( ho + Vs) <6 ° < - a)ho - - S. 

^7(1—0) V v a / v ^70(1-0) 



VTCpa) V y/fa(l-a) 

< 6»o < 00 



< /io < 



y / 7a(l-a) 

/i < 

- -r==S - h o + ±)<°o< - V7(l - a)h - , ; ■ 



(B9) 
(BIO) 
(BH) 



The minimal distances are given by 



,tz[ _ {\/ia{\ - a)6 + k + y/ah o y 

rain 



o[l + 7(1-0)] 



rf 1 ' 2 - h 2 
a min — n 



tz' 3 _ (v^yo(i - a)0 o + k- y/ah o y 

rain 



o[l + 7(1-0)] 



(B12) 
(B13) 

(B14) 



and redefining n/ y/~fa{l — a) + 8o — > 9 we find (|3l|)-(33) 



APPENDIX C: RS STABILITY 



Starting from the stability matrix formed by the sec- 
ond derivatives of 4> (recall eq.(A8)) with respect to the 
order parameters and the conjugated variables, we find 
that only transverse fluctuations are relevant. 

These transverse fluctuations are characterized by four 
eigenvalues with degeneracy n(n — 3)/2, given by the 
roots of the fourth degree characteristic polynomial P(X) 



P(A) = 



A, - A 
A c 
1 




A c 
A. r - A 

f 



[(A, - A) (A ? - A) - 1] [(A,, - A) (A? — A) — 1] — A|! ( A ? - A) ( A ? - A) (CI) 



with the coefficients A given by 



A 9 = 


a 

?J 


f V(h )V(^j6 - 


A r = 


a 

^ J 


f V(h )V{^6 a - 


A c = 


a 

qr J 


f V(h )V(^y6 - 


A ? = 
A ? = 


(1- 
(1- 


q? 

r) 2 =7 2 (l-9) 2 



r 



(C2) 
(C3) 

(C4) 

(C5) 
(C6) 



where we recall that (1 — r) = 7(1 — q) and the function 
[i\t(ho, Go) is defined in @. 
Next, the limit q — > 1 has to be taken. Using 



the asymptotic expansion of [l]{(/io,#o) discussed in ap- 
pendix B we can compute the asymptotic behavior of 
the coefficients A q , A r and A c . After a lot of algebra we 
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finally arrive at the expressions (|39|)-([l3|) with the inte- 
gration regions and minimal distances given by (|2^)-(33). 
In this limit, it turns out that an analytical expression 
can be found for the eigenvalues. First, we notice that 
the determinant of the matrix remains finite in the limit. 
Since the determinant is the product of the eigenvalues, it 
follows that this product is finite. Two possibilities arise, 
either all eigenvalues are finite, or two of them tend to 
zero and two to infinity with the same ratio. It is not 
hard to prove that the first choice is incorrect. Hence, 
two of the eigenvalues have to behave asympotically as 
(1 - q) ±n . One can check that only n = 2 is possible. 
This allows us to split P(A) into two polynomials which 



give the solutions around zero and around infinity. These 
polynomials read 



P (X) = [A,(A ? -A)-l][A r (A ? -A)-l] 
-A^(A ? -A)(A ? -A) 
-Poo (A) = (A q -X){A r -X)-A 2 c . 



(C7) 



From these two polynomials the four eigenvalues (p5|)- 
( |38| ) can be found. We remark that in the limit a — > 1 we 
find back the stability criteria for the original Gardner 
capacity problem. 
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